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Abstract 

O | We study the pseudo-Hermitian systems with general spin-coupling point 

interactions and give a systematic description of the corresponding bound- 
ary conditions for PT-symmetric systems. The corresponding integrability 
for both bosonic and fermionic many-body systems with PT-symmetric 
contact interactions is investigated. 
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^ | Self-adjoint quantum mechanical models describing a particle moving in a local singular 

potential have been extensively discussed [1, 2, 3]. The integrability of one dimensional 

many-body systems with self-adjoint contact interactions has been studied according to 

Yang-Baxter relations [4]. The results are generalized to the case of particles with spin- 
Ch ' 

c3 \ coupling interactions [5]. 

^ ! 

PT-symmetric quantum mechanical models have been studied from some mathemati- 
cal and physical considerations [6]. In [7] the classification and spectra problem of PT- 
symmetric point interactions are investigated. The integrability of many-body systems with 

hi ' 

PT-symmetric interactions is clarified [8]. The 5-type spin-coupling interactions with PT- 
symmetry is discussed in [9]. 

In this letter we study the boundary conditions for PT-symmetric point interactions of 
particles with spin-coupling, and the integrability of bosonic and fermionic many-body sys- 
tems with PT-symmetric, spin-coupling contact interactions characterized by these bound- 
ary conditions. 

One dimensional quantum mechanical models of spinless particles with point interac- 
tions at the origin can be characterized by separated or nonseparated boundary conditions 
imposed on the (scalar) wave function (p at x = 0. The family of point interactions for the 
Schrodinger operator — ^ can be described by unitary 2x2 matrices via von Neumann 
formulas for self-adjoint extensions of symmetric operators: 




where ad — bc = 1, 9, a, b,c,d G iR. ^(a;) is the scalar wave function of two spinless particles 
with relative coordinate x. (1) also describes two particles with spin s but without any spin 
coupling between the particles when they meet, in this case f represents any one of the 
components of the wave function. 

The separated boundary conditions with respect to self-adjoint interactions are described 

by 

<p'(0 + ) = h + y>(0 + ) , <p'(0_) = h-<p(0.), (2) 

where /i ± G M U {oo}. h + = oo or h~ = oo correspond to Dirichlet boundary conditions 
and h + = or h~ = correspond to Neumann boundary conditions. 

The family of PT-symmetric point interactions is described by the boundary conditions 
at the origin of one of the following two types 



with the real parameters b > 0, c > —1/b (if the parameter b is equal to zero, then the last 
inequality can be neglected), 9,<p £ [0,27r); and 

h v'(0 + ) = /i 1 e i V(0+) , W(O-) = -/ne-*V(0-) (4) 

with the real phase parameter 9 G [0, 2tt) and with the parameter h = (ho, hi) taken from 
the (real) projective space P 1 . 

For a particle with spin s, the wave function has n = 2s + 1 components. Therefore 
two particles with contact interactions have a general boundary condition described in the 
center of mass coordinate system by: 

where ip and ip' are n 2 -dimensional column vectors, A, B, C and D are n 2 x n 2 matrices. 
The boundary condition (5) includes not only the usual contact interaction between the 
particles, but also the spin couplings of the two particles. 

For self-adjoint point interactions, due to the required symmetry condition of the Schrodinger 
operator: 

d 2 d 2 

for any u, v G C°°(IR\ {0}), the matrices A, B,C, and D are subject to the following 
restrictions: 

A ] D-C ] B = E, B ] D = D ] B, A ] C = C^A, (6) 
where f stands for the conjugate and transpose. 
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Corresponding to (2) the separated boundary conditions are given by 

V/(0 + ) = G + V(0 + ), V'(O-) = G"^(0_), (7) 

where, for the self-adjoint point interactions, G ± are Hermitian matrices. 

We consider now the boundary conditions describing the point interactions with PT- 
symmetric and spin-coupling interactions. By applying PT-operation to (5), we have 

i/j \ (A* -B* W V 
i>< ) o _-{-C* D* )[ip> 

where * stands for conjugation. Hence A, B, C, D satisfy 

AA* - BC* = E, DD* - CB* = E, BD* = AB*, CA* = DC*, (8) 

where 1 is the n 2 x n 2 identity matrix. The boundary condition (5) with A, B, C, D 
satisfying (8) represents PT-symmetric point interactions with spin coupling. 

Accordingly the separated type boundary conditions with respect to PT-symmetric in- 
teractions are given by 

^(0+)=F^(0 + ), V'(0-)=Cty(0_), (9) 

where G = —F*. 

The case A = D = E,B = 0,C = C* corresponds to a Hamiltonian with PT-symmetric 
5-type interactions (when C is further symmetric, the system is both PT-symmetric and 
self-adjoint). The case A = D = H, C = 0, B = B* corresponds to a PT-symmetric 
Hamiltonian H of the form: 

H = -D 2 X {1 + B5) - BDJ', 

where D x is defined by (D x f)((p) = — f{^f), for / e C£°(lR/{0}) and if a test function 
with a possible discontinuity at the origin. H is self-adjoint when B is Hermitian [5], and 
H is both self-adjoint and PT-symmetric when B is real symmetric. 

Concerning the integrability of many-body systems with spin-coupling interactions, we 
first consider two-particle case. Let e a , a — 1, ...,n, be the basis (column) vector with the 
a-th component as 1 and the rest components 0. The wave function of the system is of the 
form 

n 

i> = H (f>a/3(xi,x 2 )e a <S> ep. (10) 

a,f3=l 

According to the statistics ip is symmetric (resp. antisymmetric) under the interchange 
of the two particles if s is an integer (resp. half integer). Let k\ and hi be the momenta 
of the two particles. In the region x\ < x%, in terms of Bethe hypothesis the wave function 
has the following form 

$ = Ul2e i(k lXl +k 2 x 2 ) + U2ie i(k 2 x 1+ k lX2 ) = Ui2e i(K 12 X-k 12 x) + U2ie i(K 12 X+k 12 x) ? ( U j 
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where X = (x\ + x 2 )/2, x = x 2 — Xi are the coordinates of the center of mass system, 
K12 = ki + k 2 , k±2 = (ki — k 2 )/2, u 12 and u 2 \ are n 2 x 1 column matrices. 
In the region X\ > x 2 , 

$ = (P l2 Ul2 )e i(Kl2X+kl2x) + (P l2 u 21 )e^ Kl2X - kl2X \ (12) 

where according to the symmetry or antisymmetry conditions, P 12 = p 12 for bosons and 
P 12 = —p 12 for fermions, p 12 being the operator on the n 2 x 1 column that interchanges the 
spins of the two particles. Substituting (11) and (12) into the boundary conditions (5), we 
get 

U12 + u 21 = A P 12 {u 12 + u 2 i) + ik 12 B P l2 {u l2 - u 21 ), 
ik 12 (u 21 - uu) = C P 12 (u 12 + U21) + ik 12 D P l2 {ui 2 - u 21 ). 
Eliminating the term P 12 u 2 i from (13) we obtain the relation 



(13) 



u 21 = Y 2 12 u 12 , (14) 

where 

Y 2 \ 2 = [(A - ikuB)- 1 - tk 12 (C - ikuD)- 1 ]- 1 

[(A - ik 12 B)-\A + ik l2 B)P 12 - (C - ik 12 D)-\C + ik 12 D)P 12 (15) 
-(A - ikuB)' 1 - ik 12 {C - ik l2 D)- 1 }. 
Similarly, with respect to the separated type boundary condition (9), we have 

ik 12 (u 21 - u 12 ) = F(u 12 + u 2 i), 
tk 12 P 12 (u 12 - U21) = -F* P 12 (u 12 + u 21 ). 
The two equations above are compatible when F satisfies 

(ik 12 - F)-\ik l2 + F) = P 12 (tk l2 - F*)-\ik l2 + F*)P 12 . (17) 

(17) is satisfied when F commutes with P 12 , which implies that F is real. For the case of 
spin- 1 (n — 2), F is generally of the form 



(16) 



F = 



I a ei e\ c \ 

e 3 f 9 e 2 

e3 9 f e 2 

\ d 64 64 b J 



(18) 



where a, b, c, d, f, g, ei, e 2 , e 3 , e 4 G M. In stead of (15), from (16) we have, upon to the 
condition (17), 

Y 2 12 = {ik l2 -F)~\ik l2 + F). (19) 

In the following we consider the integrability of systems with point interaction described 
by the boundary condition (9). For a system of iV identical particles with PT-symmetric 



contact interactions characterized by the separated type operator (19), the wave function 
in a given region, say X\ < x 2 < ... < x N , is of the form 

n 

* = <i>ai,...,a K ( x u—,x N )e ai ®...®e aN 

ai,...,ajv=l 

— u i2 N e i ( klXl+k2X2+ '" +kNXN " > + U21 N e i ( k2Xl+kiX2+ --- +k N x N) (20) 
+ (iV! - 2) other terms, 

where kj, j = 1, N, are the momentum of the j-th particle, u are n N x 1 matrices. The 
wave functions in the other regions are determined from (20) by the requirement of symmetry 
(for bosons) or antisymmetry (for fermions). Along any plane X{ = i G 1, 2, TV — 1, 
we have 

7/ = Y jj+1 v (21 ) 

u 'Oi02---OjQ!j+i---ajv aj + ia>j u, a>ia2...aj+iaj...apf! \' Lj± J 

where 

Yaj+ioj — [ik aj a j+1 — Fjj+i] 1 [ik aj a j+1 + Fjj+i]- (22) 

Here k ajCtj+1 = (k aj — k aj+1 )/2 play the role of momenta and P" +1 = for bosons and 
pjj+i — _pjj+i f or fermions, where is the operator on the n N x 1 column u that 

interchanges the spins of particles j and j + 1. Fjj +1 stands for the application of the 
operator F to the jth and j + 1th particles. 

For consistency Y must satisfy the Yang-Baxter equation with spectral parameter [10], 

ym,m+lym+l,m+2ym,m+l ym+l,m+2ym,m+lym+l,m+2 ( C )'X\ 

Y ij Y kj 1 ki — X ki 1 kj Y ij ■ K^ 6 ) 

It is straight forward to verify that the operator Y given by (22) satisfies the Yang-Baxter 
equation (23). Therefore this system is integrable with the exact wave functions given by 
(20). 

As F is a matrix satisfying (17), it's spectra are not real in general. For instance, the real 
matrix (18) gives complex eigenvalues. Hence the PT-symmetric (separated type) contact 
couplings contain non-Hermitian interactions with complex and real spectra, while the real 
spectra, e.g. spectra of real F commuting with P 12 , covers part of the real spectra from the 
Hermitian interactions, when the matrix F is further symmetric (e^ — ei, — d — c in 
the case of (18)). 

As for bound states, let us assume that F has real spectra. Let Y be the set of n 2 
eigenvalues of F. For any X a e T such that \ a < 0, there are 2 N ( N ~ 1 ^ 2 bound states for the 
iV-particle system, 

< = vae U(0(x k - x{) + e kl 6{ Xl - x k ))e Xa ^ ^\ (24) 
k>i 

where v a ^ is the spin wave function and e = {e k i : k > I}; e k i = ±, labels the 2 Ar ( Ar ~ 1 )/ 2 -fold 
degeneracy. The spin wave function v here satisfies P !3 v = e^v at for any i ^ j, that is, 
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V^Vae = tijVat for bosons and p l ^v a± = —e^v a ^ for fermions. The energy of the bound state 

E a = -^N(N 2 - 1). (25) 

We have studied the boundary conditions for systems with general PT-symmetric spin- 
coupling point interactions, and the corresponding integrability for both bosonic and fermionic 
many-body systems with separated-type PT-symmetric contact interactions (19). The scat- 
tering matrices can be also studied similar to the case of self-adjoint interaction case. The 
spectra and integrability for many-body systems with PT-symmetric contact interactions 
described by the Y operator (15) for the non-separated boundary conditions can be studied 
accordingly in terms of the Yang-Baxter equation (23), though it could be quite complicated 
as A, B, C, D subject to the conditions (8). 
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